Spin transport and magnetoresistance in F/S/F spin valves 
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We consider spin transport and spin relaxation in superconductors using the quasiclassical theory 
of superconductivity. We include spin relaxation due to spin-orbit interaction as well as magnetic 
impurities, and show that the energy dependence of the spin-flip rate is different for these two 
mechanisms. In ferromagnet-superconductor-ferromagnet systems made of Co and Al, interface 
resistances can be small compared to bulk resistances. This simplifies the description of transport in 
Co/Al/Co spin valves, for which we numerically calculate the temperature and Al length dependence 
of the magnetoresistance. 



I. INTRODUCTION 

Ferromagnetism and superconductivity are two com- 
peting phenomena in condensed matter physics. In con- 
ventional low temperature superconductors, transport of 
spins beyond the coherence length is prevented by the 
formation of spin singlet Cooper pairs. Consequently, 
due to the competing ordering of ferromagnets (F) and 
superconductors (S) in hybrid structures, many nontriv- 
ial physical effects occur 1,2 and there are interesting sug- 
gested applications such as an absolute spin valve effect 3 
and solid state memory elements. 4 

Experimental studies of F/S contacts in the diffusive 
limit showed that the resistance can both decrease and in- 
crease relative to the resistance above the critical temper- 
ature (T c ) of the superconductor. 5-7 Theoretically it was 
shown that the temperature dependence of this resistance 
depends sensitively on the contact transparency. 8-10 The 
resulting resistance is determined by an interplay be- 
tween the energy-dependent interface resistance and spin 
accumulation at the interfaces due to reduced spin trans- 
port into the superconductor. 

Transport of spins through the bulk of superconduc- 
tors was recently studied experimentally in an F/S/F 
heterostructure. 11 Here, a decreased magnetoresistance 
(MR) in the superconducting state was interpreted as a 
loss of spin memory. Theoretical work on bulk spin trans- 
port in superconductors in the inelastic regime 12 and 
the elastic regime 13 describes the reduced penetration of 
spins by spin flipping and reduced penetration of spin po- 
larized quasiparticles. The F/S/F system of Ref. 11 has 
been analyzed by assuming a spatially homogeneous su- 
perconducting order parameter and neglecting spin flip. 14 
However, a thorough understanding of spin transport in 
F/S systems requires a description of the spatially de- 
pendent order parameter in each component as well as 
the quasiparticles driven out of equilibrium. A theoret- 
ical description of the F/S/F system, where the spatial 
variation of the order parameter, energy-dependent spin 
flipping in the superconductor, and the effect of the inter- 
faces is taken into account, has to the best of our knowl- 
edge not been published. 

In order to study the bulk spin transport properties, it 



is important to have control over the influence of inter- 
faces. Typically, in spin valve structures both interface 
resistances and bulk resistivities contribute to the MR 
and are affected by superconductivity. In this paper, we 
study a superconducting spin valve system, where the 
interface resistances are negligible. In that case, a sim- 
plified treatment of the F/S boundaries is possible so that 
bulk effects can be studied independently of interface ef- 
fects. As we discuss later, a possible candidate to realize 
a spin valve with small interface resistance could be a 
Co/Al/Co system. To describe the transport through 
a superconducting spin valve, we present a theoretical 
framework that describes the spin dependent transport in 
superconductors in linear response. Spin flip scattering 
from magnetic impurities as well as spin-orbit interac- 
tion is included in our description, and the full spatial 
dependence of the pairing potential is calculated self- 
consistently. We use this formalism for numerical cal- 
culations of the magnetization-configuration dependent 
transport of a Co/Al/Co spin valve. This demonstrates 
the suppression of spin transport through the supercon- 
ductor. 

The paper is organized in the following way: Section 
II describes the equations governing elastic transport in 
a diffusive superconductor. Section III outlines the spe- 
cific geometry studied and the approximations used. In 
section IV we discuss the numerical results. Section V 
summarizes and concludes our work. 



II. TRANSPORT THEORY 

Using the Keldysh theory in the quasiclassical approx- 
imation, we have previously derived kinetic equations for 
transport of charge, energy, spin and spin energy in diffu- 
sive, s-wave superconductors in the presence of spin-flip 
scattering by magnetic impurities. 13 We will now sup- 
plement that treatment with expressions for spin-orbit 
induced spin relaxation, and derive the resulting trans- 
port equations in linear response. For an explanation of 
the notations used below we refer to Ref. 13. 
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The spin-orbit interaction Hamiltonian is 

a' a J 

(1) 

where 7 is the interaction strength, r is the vector of 
Pauli matrices and V5 mp is the impurity scattering po- 
tential. The spin-orbit contribution to the self-energy in 
the Eilcnbcrgcr equation 15 is 
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where we have defined the spin-orbit scattering time 
1/tso = 8j 2 pp/9r. Here pf is the Fermi momentum, r is 
the elastic scattering time, a is a vector of 4x4 matrices 
with the Pauli matrix and its transpose on the diago- 
nal block, i.e. a = diag(r, f T ), ps = diag(l, 1, —1, —1), 
and g s is the isotropic part of the Green's function in 
Keldysh-Nambu-spin-space. Using a convenient repre- 
sentation of the Green's functions, we obtain equations 
that determine the distribution functions and currents. 

Linearized kinetic equations for charge transport in 
diffusive superconductors were obtained by Schmid and 
Schon 16 and have been successfully applied to describe 
various transport phenomena. To study spin-dependent 
transport it is necessary to include equations that de- 
termine the spin-current. The relevant equations in the 
linear response regime are developed below. The approx- 
imations are valid when deviations from equilibrium val- 
ues are small. We also assume that any static super- 
current is small, i.e. that there is no Joscphson effect. 
The transport theory is formulated in terms of the phys- 
ical particle and energy currents (including particles and 
holes). These are given by the distribution functions /it 
and /i L 17 and the spin resolved functions /its and /ils> 
as well as generalized diffusion coefficients Dt, £>l and 
renormalization factors «tt, Q^tsts f° r relaxation pro- 
cesses. The spin-resolved distribution functions can be 
expressed by the particle distribution function as 

; T (£)-; X (£) ; T (-£)-;x(-g) 
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The spectral (retarded) properties depend on the com- 
plex function 9(E) which is determined by the so-called 
Usadel equation. 17 To describe spin polarized transport 
in voltage biased systems in linear response, it unnec- 
essary to calculate /il and /ilSj so the equations that 
determine these functions have been omitted below. 

The charge current and the spin current in S are given 
by integrals over the spectral quantities. The charge and 
spin current carried by quasiparticles is 29 

1 f°° r)h 
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where A is the area of the wire and Nq is the density 
of states at the Fermi level for both spins in the normal 
state. Additionally, charge current is carried by the su- 
percurrent, so that the total charge current is constant. 
The distribution functions /it and /its are determined 
by the diffusion equations 
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Here r m is the spin flip scattering time due to magnetic 
impurities and r so the spin flip scattering time due to 
spin-orbit coupling, both evaluated in the normal state. 
The spectral quantities are given in terms of 9(E, x). Wc 
compute that the renormalization of the scattering rates 
are 
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The effect of spin-flip scattering by spin-orbit interaction 
with renormalization factor a^grpg is a new result that 
did not appear in our previous article. 13 Its renormaliza- 
tion is different from the renormalization of the spin-flip 
scattering by magnetic impurities. The complex function 
9 is determined by the Usadel equation 
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Note that the spin flip term in (9) arises from magnetic 
impurities only since spin-orbit scattering does not lead 
to pair breaking and consequently docs not influence the 
spectral properties of the superconductor. This equation 
must be solved in conjunction with the self-consistency 
relation 
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where A is the electron-electron interaction strength and 
Ed the Debye cut-off energy. 

An applied voltage is taken into account as a bound- 
ary condition for the distribution functions, /it and /its- 
In a reservoir with electrochemical potential \i we have 
in linear response the equilibrium distributions h\ = 
-(3/j,/ (2cosh 2 (/3£;/2)) and h^ s = 0. 

The different renormalization factors c^tsts an d a TSTS 
arise from spin flipping by spin-orbit interaction or mag- 
netic impurities. In general, a^g-jg and aJpgrpg depend on 
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the spectral properties of the superconductor through 8. 
In the BCS limit, valid for large bulk superconductors, 
the energy dependence of these factors is completely dif- 
ferent and correspond to the so-called type-I or type-II 
coherence factors. 18 ' 19 Using the BCS solution for the 
Green's functions we find that for energies below the 
gap (for which there are no quasiparticles in the BCS 
limit), both cc^sts anc ^ a TSTS van i sn 7 an d above the gap 
a s ° STS = 1 and a£ STS = (E 2 + A 2 )/(E 2 - A 2 ) > 1. Fur- 
thermore, we see from Equations (8b) and (8c) that for 
any 8(E), «xsts < a TSTS- This i m P nes that for a given 
normal state spin-flip length, the rate of spin flipping in 
the superconducting state is higher when the dominant 
spin-flip scattering mechanism is caused by magnetic im- 
purities than if it is caused by spin-orbit interaction. 

III. MODEL 

We consider an F/N/S/N/F wire as shown in Figure 1. 
It is assumed that the ferromagnets (Co) are connected 
via normal metals (Cu) to the superconductor (Al). The 
distribution functions in the ferromagnet and the nor- 
mal metal are determined by the Valet-Fert transport 
theory 20 and in the superconductor by the theory de- 
scribed in the previous section. An applied bias causes 
spin polarized quasiparticles to be injected into the S 
layer. We assume that the magnetization of the F parts 
are either parallel (P) or antiparallcl (AP). Because of 
renormalized spin flip rates and a reduction of the gen- 
eralized spin diffusion coefficient (-Dl) in the supercon- 
ductor, the magnetoresistance MR = (i? AP — R p )/R p 
is reduced for temperatures below T c compared to the 
normal metal state. 




FIG. 1: Diffusive wire consisting of F,N and S elements. 

In order to determine the dominant contributions to 
the resistance of the system, we examine the magnitude 
of the resistance of the F/N interface (-R f /n) compared 
to the bulk resistance in F within a spin- flip length (R F ). 
The latter quantity is the largest resistance of the ferro- 
magnet within its spin-active part. To this end, consider 
the ratio 

Rf/n _ AR F / N 
We assume that F layers are made of Co, N layers of 



Cu and the S layer of Al. The bulk resistance of Cu 
as well as the proximity effect is neglected since the Cu 
layer is very thin, and in addition the typical resistiv- 
ity of Cu is smaller than that of Co or Al. The inter- 
face resistance (AR F / N ), resistivity (p F ) and spin diffu- 
sion length (iff) for Co is reviewed in Rcf. 21. It is 
found that ^4i?cu/Co ~ 0.5 ffim 2 at 4.2 K where A is 
the cross section area. The renormalized resistivity 2 is 
p* Co ~ 75 nftm, and Z s c f ° = 59 nm at 77 K. 22 Thus we 
can conclude that for Co R F / N /R F { ~ 0.1 as a least es- 
timate since the spin diffusion length should be longer 
at 4.2 K. This means that it is a valid approximation to 
disregard the interface resistances for long enough sam- 
ples. The N/S interface resistance between Cu and Nb 
above the critical temperature is found to be larger than 
the F/N resistance, 23 ^4i?cu/Nb ~ 1-10 ffim 2 , and would 
give i?N/s/-Rff ~ 0.2. With Al as the superconducting 
layer we expect no higher interface resistance. We may 
also argue that the bulk resistance for dirty Cu/Co lay- 
ers scales as AR^° lk ~ 0.1L[nm] film 2 , 24 where L is the 
length of the layers expressed in nm. Thus the bulk re- 
sistance for a slice of length l s f should be much larger 
than the interface resistance. Note that in direct F/S in- 
terfaces the reported interface resistance is considerably 
higher AR Nb/Co ~ 3 ffim 2 . 23 

The estimates above show that the interface resistances 
are much smaller than the relevant bulk resistances with 
the materials chosen here. We will later check that 
the change in resistance from normal to superconducting 
state is larger than the interface resistances. A possible 
approximation is therefore to neglect the interface resis- 
tances. This allows us to effectively do calculations for 
an F/S/F system with the boundary condition that the 
generalized diffusive current should be continuous which 
implies that the function 8 is continuous at the interface. 
For strong ferromagnets the superconducting proximity 
effect into the ferromagnet is negligible and we have by 
continuity 8 — > in the superconductor close to the F/S 
interfaces. This means e.g. that the gap vanishes at 
the interface. In this case, it is the bulk transport prop- 
erties that dominate the system, and there are no free 
parameters so that it is possible to give an unambiguous 
description of the transport properties. This is our aim 
in the rest of the paper. 

The F/S/F system as shown in Fig. 1 was motivated 
by the experiments of Gu et al. However, in those 
experiments Py was used for the ferromagnet, and be- 
cause of the very short spin diffusion length in this al- 
loy (Z^ y = 5.5 nm) the interface resistances are of the 
same order as R F f . Consequently, in these experiments 
both the interface resistances and the bulk resistance of 
Al are governed by superconductivity. Thus the model 
discussed above is not applicable, and the resistance of 
the spin polarizing interface must be taken into account. 
To be specific, we no longer have that 8 — > at the 
interfaces, and superconductivity is not completely sup- 
pressed at the interface as in the Co/Al/Co system. Us- 
ing the approximations discussed above in calculations 
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for the Py/Nb/Py system of Ref. 11 would therefore 
give a too low T c . Numerical simulations and compari- 
son with Ref. 11 show that this is indeed the case (not 
shown). A complete description of this experiment re- 
quires boundary conditions for the spin polarizing inter- 
faces given by scattering theory. This would describe the 
proximity effect in N as well as a reduction of the su- 
perconducting pairing amplitude close to the interface. 
However, as noted by Huertas-Hernando et al, this ap- 
proach would require full knowledge about the interface 
scattering matrix, which is generally not available except 
for in simplified models at this moment. 

IV. CALCULATIONS 

We have performed numerical calculations for a 
Co/Cu/Al/Cu/Co spin valve. Parameters for the super- 
conductor are mostly taken from Ref. 26. The bulk 
value of the pairing potential at zero temperature is 
Ao = 192 yueV and the critical temperature T c = 1.26 K 
with interaction parameter NoX/2 — 0.18. 27 The normal 
state diffusion coefficient of Al is D = 160 cm 2 /s, and the 
density of states at Fermi level N = 2.2 • 10 47 J^m" 3 
corresponding to a resistivity p^i = 11 nfhn. The nor- 
mal state spin-flip relaxation length by spin-orbit inter- 
action is given by the sample independent parameter 
e = l s {/l i=a 30, 28 and we assume that the elastic mean free 
path is I — 37 nm. This gives l s { = 1.11 /um for spin-orbit 
induced spin-flip. In calculations for magnetic impurity 
induced spin-flip we take the normal state value of the 
spin-flip length identical to the spin-orbit induced l s {, but 
in general this length is determined by the impurity con- 
centration which is sample specific. We take the length of 
the (identical) ferromagnetic elements to be 100 nm with 
a bulk spin asymmetry f3 = OA. 21 Figure 2 shows the 
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FIG. 2: Spatial variation of the pairing potential at T/T c = 
0.40 scaled by Ao. Inset: Density of states at positions 0.0, 
0.2 and 0.6 fim into the S wire. The curve evaluated at x = 
0.0 fj,m is identical to the normal state DOS (flat curve). 

spatial variation of the pairing potential resulting from 



complete suppression of superconductivity at the F/S in- 
terfaces at reduced temperature T/T c = 0.40 for a 1.2 /im 
Al wire with magnetic impurities. The density of states 
at various locations in the superconductor is shown in 
the inset, and resembles the bulk BCS shape close to the 
center of the wire where the gap is largest. 
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FIG. 3: Temperature dependence of resistance in parallel 
geometry for spin-orbit (solid lines) and magnetic impurity 
(dashed lines) induced spin flip. The thick curves are with Al 
length 1200 nm, and the thin curves with Al length 900 nm. 
T c is lowered by the presence of the magnetic impurities. 

A calculation of the resistance of the F/S/F system 
for parallel magnetizations is shown in Fig. 3. The AR P 
values above the critical temperature agrees with analyt- 
ical results based on the Valet-Fert theory. Below T c the 
resistance drops rapidly, but remains finite in the limit 
T — > 0. The change in resistance from normal to super- 
conducting state is of the order of 2 - 6 film 2 depending 
on the length of the superconductor, and this change is 
larger than the typical interface resistance, which should 
be checked as noted in Section III. The resistance of the 
system below T c is due to the F elements as well as the 
regions in the S wire next to the F/S interface where 
there is conversion of current into supercurrent. 26 The 
systems with magnetic impurities has the higher resis- 
tance as T — ► since the length of the resistive region 
near the interfaces is longer. This is because the conver- 
sion of current into supercurrent happens over a length 
scale determined by the coherence length £ = ^/HD/2irA 
which for a superconductor with magnetic impurities is 
longer since A is suppressed due to a term in the Usadel 
equation (9). 

The dependence of the resistance on the magnetiza- 
tion configuration is shown in Fig. 4 where the excess 
resistance AR = A(R AP — i? p ) is plotted as function 
of temperature. We show curves for systems with only 
spin flip scattering from magnetic impurities or spin-orbit 
interaction. The systems with magnetic impurities pro- 
vide a weaker suppression of the spin signal than sys- 
tems with spin-orbit interaction. The opposite could be 
expected since as noted above a™ STS > a™ STg . How- 
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FIG. 4: Temperature dependence of the magnetoresistance for 
spin-orbit (solid lines) and magnetic impurity (dashed lines) 
induced spin flip. The top set of curves is a system with Al 
length 900 nm, and the lower curves with Al length 1200 nm. 



ever, the pairing potential is lower in a superconductor 
with magnetic impurities due to the detrimental effect 
of the impurities on superconductivity, and this is the 
dominant effect. This is confirmed by simulations of sys- 
tems with equal strengths of the pairing potential, in 
which magnetic impurities is the strongest spin relax- 
ation mechanism. From Fig. 4 we see that the difference 
in suppression of spin signal between spin-orbit and mag- 
netic impurity induced spin-flip is smaller for the longer 
wires, since in this case the difference in A is also smaller. 
For long wires the excess resistance tends to zero at low 
temperatures as expected, because in this case the trans- 
port of spins through the superconductor is completely 
suppressed. 

In Fig. 5 we show the spatial variation of the quasi- 
particle charge and spin current and supercurrent for the 
F/N/S/N/F spin valve with parallel magnetizations. The 
charge current is constant in the F parts of the wire, and 
is gradually converted into supercurrent in S. Spin cur- 
rent injection into S is suppressed, as a comparison with 
the magnitude of spin current in the normal state shows. 
This leads to spin- accumulation in F at the interfaces. 
We see that the spin current is reduced below T c inside 
the superconductor due to Cooper pairing. On the other 
hand, the total charge current increases below T c due to 
the reduced resistance of the superconductor. In Fig. 6 
we show the spin accumulation /if — jj,^ for the same sys- 
tem. Comparison with the normal state shows that the 
spin accumulation is larger in the S case, due to the re- 
duced penetration of spins into S and since the net spin 
current out of the reservoirs is larger in the S case because 
the total resistance is lower. The spin accumulations that 
build up at the interfaces are relaxed through spin flip in 
S. These spin accumulations can be measured e.g. by 
tunnel coupling between the superconductor and a third 
probe ferromagnet. 

Qualitatively, our results for the MR are in agreement 
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FIG. 5: Top panel: Spatial dependence of the spin current for 
the F/S/F structure with Al length 1.2 fj,m at T/T c = 0.40 
(spin-orbit induced spin- flip). The normal state (T > T c ) 
spin-current is shown with dotted line. The F/S interfaces 
are at x — /im and x = 1.2 /im. Bottom panel: Spatial 
dependence of the quasiparticle charge current (solid line) 
and supercurrent (dashed line) for the same system as in top 
panel. Normal state current is shown with dotted line. 
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FIG. 6: The spatial dependence of the spin potential for the 
F/S/F structure with Al length 1.2 fim at T/T c = 0.40 (solid 
line) and T > T c (dotted line) (spin-orbit induced spin-flip). 
The F/S interfaces are at x = /j,m and x = 1.2 pirn. 



with the experiment by Gu et al. A contribution from 
the interfaces which is most probably important in the 
experiment, will not qualitatively change the properties 
of the system except for a higher T c as noted above. 
Therefore, quantitative differences between the experi- 
ment and our calculations using material parameters cor- 
responding to the system in Ref. 11 are not surprising. 
A more detailed theoretical analysis, which accounts for 
interface resistance, should be made to enable a quantita- 
tive comparison with the experiments of Ref. 11, but this 
is beyond the scope of our present work. We emphasize 
again, that our predictions are, however, experimentally 
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testable in Co/Al/Co spin valves, which can be fabricated 
using state-of-the-art technology. 

V. SUMMARY AND CONCLUSIONS 

In conclusion we have studied spin transport properties 
of an F/S/F trilayer. We have developed transport equa- 
tions using the quasiclassical theory of superconductivity 
and included the effects of spin-flip relaxation. An ex- 
perimental system is proposed where interface resistance 
can be neglected and a simple description of the physics 
is possible. For this system we have performed numer- 
ical calculations of the magnetization-configuration de- 



pendent resistance. This demonstrates the dependence 
of the spin transport suppression on different spin flip 
mechanisms, i.e. magnetic impurities and spin-orbit in- 
teraction. 
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